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ABSTRACT.  This  report  describes  the  results  of  testing  a  series  of 
composite  and  double-base  propellants  in  a  T-burner  facility.  The 
investigation  reveals  the  unstable  combustion  behavior  that  is  suscep¬ 
tible  of  excitation  under  acoustic  pressure-coupled  conditions.  The 
degree  of  instability  is  quantitatively  described  in  terms  of  the  re¬ 
sponse  function,  which  is  a  well-defined  characteristic  of  the  propel¬ 
lant.  The  effect  of  pressure  and  frequency  has  been  systematically 
investigated  for  the  frequency  range  of  500-10,000  cps  and  for  pres¬ 
sures  up  to  nearly  3,000  psia.  Such  an  extended  coverage  has  permitted 
the  determination  of  a  very  general  pattern  of  behavior,  including  the 
maximum  instability  conditions,  the  isoresponse  function  contours,  and 
the  neutral  boundary  separating  the  unstable  and  stable  combustion 
regimes  in  the  1.5-inch  T-burner. 
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NOMENCLATURE 


Velocity  of  sound 

Acoustic  energy  in  a  column  of  unit  cross-sectional  area  and 
length  L 

Rate  of  change  of  E  with  respect  to  time 
E  during  growing  oscillations  with  attenuation 
Acoustic  energy  per  unit  volume 
e  during  decaying  oscillations 
Frequency 

Length  of  T-burner 

Mass  flow  rate  per  unit  area 

Amplitude  of  pressure  oscillations  at  time  t  «  0 

Amplitude  at  time  t  =  0  for  case  of  oscillations  growing  in 
presence  of  attenuation 

Linear  regression  rate  of  the  burning  surface 
Time 

Mean  flow  velocity  just  outside  the  combustion  zone 

Real  part  of  the  acoustic  admittance  (Eq.  16) 

Effective  admittance  (real  part)  including  contribution  due  to 
convection  through  the  burning  surface 

Y/ (U/P) 

Y ' / (U/P) 

Exponential  growth  rate  of  oscillations  in  absence  of  attenuation 


Y  Ratio  of  specific  heats  of  product  gas 
AP  Amplitude  of  pressure  oscillations 

AP^  Amplitude  of  pressure  oscillations  during  period  of  growth 

AP'  Amplitude  at  specific  time 

AP"  Amplitude  at  one  cycle  after  time  for  which  AP  =  AP 1 
AU  (AU)m  cos  0 

(AU)^  Amplitude  of  velocity  oscillations  just  outside  the  combustion 
zone  (velocity  perpendicular  to  the  burning  surface) . 

e  Amplitude  of  pressure  oscillation  divided  by  mean  pressure 

0  Phase  of  the  velocity  oscillation  relative  to  the  pressure 

oscillation 

V  Amplitude  of  in  phase  component  of  mass  rate  oscillation  divided 
by  mean  mass  rate 

p  Density  of  combustion  product  gas 

p  Density  of  the  solid  propellant 
s 

a  Amplitude  of  in  phase  component  of  density  oscillation  divided 
by  mean  density 

SUBSCRIPTS 

1  Pertaining  to  period  of  growing  oscillations 

2  Pertaining  to  period  of  decaying  oscillations 

m  Pertaining  to  period  of  constant  oscillations 
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INTRODUCTION 


The  Problem 

The  combustion  of  solid  propellants  in  rocket  motors  is  often 
accompanied  by  a  phenomenon  known  as  unstable  burning,  or  combustion 
instability.  This  phenomenon  can  manifest  itself  in  several  different 
ways,  depending  upon  the  type  of  propellant,  the  configuration  of  the 
motor,  and  the  operational  pressure  (Ref.  1).  Most  cases  of  combustion 
instability  result  in  motor-pressure  fluctuations  which  have  a  detri¬ 
mental  effect  on  motor  performance.  The  advent  of  very  large  motors, 
together  with  the  development  of  very-f ast-buming  propellants,  made 
combustion  instability  an  even  more  critical  problem.  Some  reliable 
criteria  for  rating  combustion  instability  are  therefore  needed. 

Among  the  types  of  combustion  instability  most  frequently  impairing 
motor  performance,  that  defined  as  "acoustic  instability"  is,  as  a  rule, 
the  hardest  to  eradicate.  It  is  encountered  when  the  motor  pressure 
oscillates  at  a  frequency  equal  to  or  near  that  of  any  of  the  acoustic 
standing  waves  of  the  motor  cavity,  either  in  longitudinal  or  trans¬ 
verse  modes.  The  peak-to-peak  amplitude  of  the  acoustic  pressure  oscil¬ 
lations  that  develop  in  very  severe  cases  of  instability  can  reach 
levels  comparable  to  the  mean  motor  pressure.  The  engineering  approach 
to  control  of  such  behavior  has  mainly  consisted  of  the  application  of 
mechanical  methods  to  reduce  the  amplitude  of  the  acoustic  pressure 
oscillations  to  tolerable  limits,  disregarding  the  basic  causative  prin¬ 
ciples  of  the  problem. 

The  dynamic  interaction  of  the  acoustic  field  with  the  propellant 
combustion  is  twofold,  due  to  the  pressure  and  velocity  components 
associated  with  the  sound  field.  The  nature  of  this  mutual  interaction 
is  highly  involved  because  of  the  peculiarities  of  the  acoustic  field 
prevailing  in  the  cavity  and  because  of  the  complexity  of  the  physical 
and  chemical  processes  occurring  during  the  combustion  of  the  propellant. 
In  order  to  make  the  problem  tractable,  the  nature  of  the  combustion  and 
sound-field  interaction  had  to  be  resolved  into  its  basic  elements,  and 
thus  the  approach  was  taken  of  studying  the  "pressure-coupled"  and 
"velocity-coupled"  contributions  separately.  Practical  experience, 
however,  has  suggested  that  the  contribution  through  pressure  coupling 
is  by  far  the  stronger  of  the  two  in  the  frequency  range  above  500  cps 
(Ref.  1). 
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Much  effort  has  been  devoted  to  the  development  of  burner  systems 
in  which  either  attribute  of  acoustic  combustion  instability  could  be 
studied  independently  of  the  other.  Considerable  success  has  been 
achieved  in  developing  a  burner  system  that  exhibits  pressure-coupled 
instability  spontaneously,  uncomplicated  by  the  presence  of  velocity 
coupling.  The  resulting  laboratory  design  is  known  as  the  T-burner, 
and  the  technique  of  testing  has  received  extensive  study  (Ref.  2). 

The  Investigation 

The  main  objective  of  the  investigation  was  to  experimentally 
determine  the  acoustic  pressure-coupled  unstable  combustion  behavior  of 
a  selected  group  of  composite  and  double-base  propellants,  using  the 
existing  T-burner  facility  at  the  Naval  Weapons  Center  (NWC) .  Emphasis 
was  placed  on  the  effect  of  pressure  and  frequency  to  try  to  detect  as 
much  of  an  overall  mode  of  behavior  as  possible.  Experimentation  at 
increasingly  higher  pressures  (up  to  2,400  psia)  over  the  available 
frequency  range  (500  to  10,000  cps)  has  led  to  the  establishment  of  a 
very  general  pattern  of  behavior  which  had  so  far  remained  unknown. 
Propellants  that  exhibit  instability  in  rocket  motors  can  generally 
be  characterized  in  the  1.5  inch  diameter  T-burner  over  the  ranges  of 
frequency  and  pressure  of  importance.  An  exception  is  the  class  of 
aluminized  propellants,  which  propellants  are  less  suited  for  study 
because  they  usually  become  stable  at  frequencies  in  the  effective  test 
range  of  this  burner  (above  500  cps)  due  to  the  stabilizing  effect  of 
AI2O3  smoke. 

The  principal  information  obtained  from  the  T-burner  tests  is 

1.  The  susceptibility  of  the  propellant  to  pressure  coupled 
instability,  as  a  function  of  frequency  and  mean  pressure 

2.  The  stability  limits  of  the  T-burner 

Various  propellant  compositions  have  been  studied  in  order  to 
compare  their  behavior.  Only  a  few  typical  double-base  propellants 
were  tested  (N-4,  N-5,  JPN,  and  X-14)  .  The  main  objective  has  been  to 
determine  the  effect  of  burning  rate  through  the  linear,  mesa,  and 
plateau  modalities  and  to  establish  the  major  differences  of  behavior 
between  double-base  and  composite  propellants.  The  objective  with  the 
composite  propellants  has  been  primarily  to  determine  the  effect  of 
(1)  ammonium  perchlorate  particle  size  and  (2)  combustion  catalysts. 

A  secondary  objective  has  been  to  determine  the  effect  of  (1)  type  of 
binder,  (2)  type  of  oxidizer,  and  (3)  amount  of  aluminum.  Some  of  the 
compositional  influences  appear  to  be  strong  and  individualistic,  and 
although  certain  trends  are  clearly  discernible,  a  formulistic  categor¬ 
ization  is  not  possible  at  present.  This  report  does  not  contain  the¬ 
oretical  work  or  a  review  of  theories;  the  brief  notes  in  "Fundamentals 
of  the  T-Burner  Method"  are  for  conceptual  definition  only. 
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FUNDAMENTALS  OF  THE  T-BURNER  METHOD 
Character  of  the  Experiment 

The  unstable  combustion  behavior  of  solid  propellants  under  acous¬ 
tic  pressure-coupled  conditions  was  investigated  in  a  T-burner.  The 
standard  design  of  this  laboratory  burner  has  a  double-end-burning  con¬ 
figuration  and  is  vented  at  the  center  to  a  surge  tank  to  afford 
constant-pressure  control  during  the  test  (Ref.  4).  The  location  of 
the  propellant  specimens  and  the  axial  mode  of  oscillation  of  the  burner 
ensure  that  the  coupling  between  the  combustion  zone  and  the  acoustic 
field  is  pressure-coupled  only.  In  the  standard  design,  the  propellant 
specimens  have  the  same  cross-sectional  area  as  the  burner.  Thus,  the 
burner  exhibits  one-dimensional,  self-excited  pressure  oscillations 
which  are  interpreted  in  terms  of  linear  acoustics  and  treated  under  the 
stipulation  that  all  the  interaction  between  the  combustion  zone  and  the 
acoustic  field  occurs  via  the  first  longitudinal  acoustic  standing-wave 
mode  of  the  burner  cavity.  A  collective  review  on  the  theory  and  use  of 
the  T-burner  technique  can  be  found  in  Ref.  5. 

In  general,  the  intensity  of  acoustic  pressure-coupled  instability 
will  depend  upon  the  combined  effect  of  several  variables: 

1.  The  amount  of  energy  released  at  sites  in  the  motor  cavity  that 
could  become  antinodes  for  acoustic  pressure  oscillations 

2.  The  modes  of  energy  released  at  the  combustion  zone  through 
kinetic  and  diffusion  processes  and  the  corresponding  reaction  times 
associated  with  these  processes 

3.  The  acoustic  damping  or  attenuation  properties  of  the  propel¬ 
lant  and  combustion  products 

4.  The  acoustic  damping  or  attenuation  properties  of  the  motor 
cavity. 

The  type  of  propellant  and  the  pressure,  frequency,  and  motor  geometry 
will  determine  in  each  case  the  magnitude  of  the  above  variables  and 
the  degree  of  their  interaction.  T-burner  experiments  cannot  resolve 
the  individual  contribution  of  each  variable,  since  they  essentially 
show  the  "global"  amplification  and  attenuation  as  a  function  of  pres¬ 
sure  and  frequency.  To  properly  elucidate  the  amplification  mechanism, 
which  is  mainly  dependent  upon  amount  and  mode  of  energy  release,  and 
the  attenuation  mechanism,  which  is  mainly  dependent  upon  particulate 
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damping  by  the  combustion  products  and  cavity  damping,  techniques  of 
greater  versatility  than  the  present  one  will  have  to  be  applied. 

Propellant  Combustion  Amplification 


The  purpose  of  this  experiment  was  to  determine  (1)  the  propellant 

which^t"  amplificatio?  cor>stant  a  and  (2)  the  limiting  conditions  at 
tion*  S  Uu  cancels  amplification.  Both  were  determined  as  fanc¬ 
ies"*  °  heC°:  r  r  pressure  and  oscillation  frequency.  In  unstable 

exponent iallv°frn^C  oscillations  appear  spontaneously  and  grow 

exponentially  from  a  small  pressure  disturbance.  If  a  system  with  no 

pressure  oscniatr1"^  ldeal  case)  *  the  exponential  growth  of  the 
pressure  oscillations  would  continue  indefinitely  until  burnout  and 

case  thramnUrS^'f^r63^6"’  35  lndlcated  in  Fi8*  K«>.  In  such  a 
theamPlitude  of  the  pressure  oscillations  AP,  as  a  function  of 
time,  would  be  given  during  the  period  of  growing  oscillations  by 

AP  *  Pe  (1) 
andaisSydLeig:aitedabTnUatl0n’  **  ^  ^  °f  °SCillati-  *  1-. 


whl r h 3^p^h 6 * ^ the  amPHtude  of  oscillations  is  limited  by  nonlinearities 

period  (Fie  lb f°nStant-“nplltude  oscillations  after  the  initial  growth’ 
P  riod  (Fig.  lb).  Upon  burnout  of  the  propellant,  the  oscillations 
decay  in  an  approximately  exponential  manner 

An  “2*2 
ap2  =  P2e 

In  interpreting  the  test  results,  it  is  assumed  that  the  attenuation 
effects  responsible  for  o2  are  the  same  as  the  atte„natlonprese“tdur- 

in  ”hich  —  the --- 

a  =  ai  +  ao 


1  This  assumption  has  not  been  systematically  evaluated  and  is 
d££„s  ^  “  S°Te  °£  «ror,  particularly  uid«  test  eon- 

ll»lts) .  I^partlcular^the^timewis^varlation 

phase  products  a" 


AP“/AP'  =  eo/f 


(a)  System  with  no  acoustic  losses  (a2  =  0) 

AP( «  p‘e«i» _  |  AP2  *  (AP)m®"a2 

w>m 


(b)  System  with  acoustic  losses  (a  >  a 2) 


Stable  test 

(c)  System  with  acoustic  losses  (a <ct£) 

a,  «  0 


FIG.  1.  Schematic  Representation  of  Oscillatory 
Behavior  in  T-Burner. 
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Figure  1(c)  shows  the  situation  when  attenuation  overrides  amplification; 
the  outcome  is  stable  burning.  The  pressure  and  frequency  conditions, 
for  a  given  T-burner,  at  which  o  exactly  matches  02  (i.e.,  04  becomes 
zero)  define  the  neutral  boundary  separating  the  stable  and  unstable 
regimes;  this  boundary  is  a  property  of  the  system  and  not  of  the  pro¬ 
pellant  alone.  Conversely,  it  should  be  emphasized  that  o  is  an  in¬ 
herent  property  of  the  propellant.  It  follows  from  Eq.  1,  as  illustrated 
in  Fig.  1(a),  that  a  can  also  be  expressed  in  terms  of  the  pressure 
amplitudes  over  1  cycle  of  oscillation: 


o  -  f  [In  (AP'VAP')] 


(3) 


where  f  is  the  frequency  of  the  pressure  oscillations.  The  word 
"constant"  used  in  connection  with  a,  04,  and  02  is  somewhat  misleading 
and  should  be  taken  with  reservation.  For  a  given  burner  and  propellant, 
these  quantities  are  assumed  to  be  constant  only  at  a  fixed  pressure  and 
frequency  and  when  considered  in  the  context  of  the  defining  relation¬ 
ships;  it  will  be  later  seen  that  they  have  a  strong  and  definite  para¬ 
metric  relationship  to  the  related  variables. 

Acoustic  Energy  Levels 

Unstable  burning  conditions  imply  the  appearance  of  certain  acoustic 
energy  levels  and  a  direct  relationship  exists  between  these  levels  and 
a.  Referring  to  the  T-burner  experiment,  as  illustrated  in  Fig.  1(b), 
the  growth  of  the  amplitude  of  the  pressure  oscillations  APj ,  as  a  func¬ 
tion  of  time,  is  given  by 


APx  =  p'eait  (4) 

where  p'  is  the  initial  pressure  disturbance  in  the  burner  cavity,  from 
which  the  spontaneous  development  of  the  full-scale  pressure  oscillations 
is  accounted  to  originate  (with  t  ■  0  at  the  time  of  the  disturbance  p'). 
The  nature  of  this  pressure  disturbance  is  left  unidentified,  although 
the  level  of  turbulence  in  the  burner  cavity,  once  propellant  burning  is 
under  way,  might  adequately  justify  the  disturbance.  In  any  case,  the 
presence  of  this  pressure  disturbance  does  not  necessarily  have  prior 
connection  with  the  oscillatory  behavior  that  it  helps  to  trigger. 

APi  attains  a  maximum  value  (AP)m  when  nonlinearities  develop  and  a 
balance  is  reached  with  no  net  change  in  acoustic  energy  after  time 
t  >  tm,  or  roughly 


(AP)  *  p'e01!1^ 
m 


(5) 


After  burnout,  the  decay  of  the  amplitude  of  the  pressure  oscillations 
AP2,  as  a  function  of  time,  follows  the  relation 


AP2  -  (AP)  e““2t 
m 


(6) 


For  a  standing  sound  wave  in  a  closed  cavity,  the  acoustic  energy  e  per 
unit  volume,  at  isentropic  conditions,  is  given  by 


e 


p(AU)2 

—4 


(AP) 2 
_ s 

4pc2 


(7) 


where  p  is  the  gas  density,  (AU)S  is  the  antinodal  amplitude  of  the  gas 
velocity  oscillations,  (AP)S  is  the  antinodal  amplitude  of  the  gas 
pressure  oscillations,  and  c  is  the  speed  of  sound.  Considering  a 
T-burner  of  length  L,  the  acoustic  energy  per  unit  cross  section  area 
during  the  growth  of  the  pressure  oscillations  will  be  denoted  by 


(AP!)2L 

ejL  *  Ej  ■  - 5 —  (8) 

4pc 


If  aj  in  Eq.  4  is  assumed  to  remain  invariant  during  the  growth  of  the 
pressure  oscillations,  Eq.  4  may  be  used  with  Eq.  8  to  write  the  specific 
input  rate  of  acoustic  energy  in  the  system  as 


1_ 

El 


AE 

At" 


El 

E!  "  2tXl 


(9) 


Under  similar  assumptions,  the  corresponding  specific  output  rate  of 
acoustic  energy  in  the  system  during  the  decay  of  the  pressure  oscil¬ 
lations  will  be 


(ap2)2l 

e2L  =  E2  -  “2—7 


and 


1_  AEz 
E2  At 


E2 

e2 


2a  2 


(10) 


(11) 


therefore,  the  specific  input  rate  of  acoustic  energy  in  the  system  due 
to  combustion  will  become 
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1  1  E 

17  “  17  =  E  =  2  (°1  +  “2)  "  2“  (12) 

or  twice  the  value  of  a,  i.e.,  a  can  be  used  as  a  measure  of  the  ampli- 
fication  of  pressure  waves  by  the  combustion.  Experimentally,  a  would 
be  obtained  from  observed  values  of  04  and  a2,  in  combination  with 
Eq.  12. 

Acoustic  Admittance 

In  most  rocket  motor  situations,  the  combustion  processes  are 
localized  near  the  propellant  surface.  Under  this  condition,  the  com¬ 
bustion  layer  can  be  viewed  for  the  purposes  of  stability  analyses  as  an 
izpHSh  ally  aCtlVe  Surface*  whose  acoustic  properties  can  be  character- 

Z  rJuan^C°UT  admittance*  The  value  of  the  admittance  depends  on 
the  combustion  dynamics,  which  in  principle  can  be  determined  either 
rom  eMlyticai  models  of  the  combustion  dynamics  or  from  experiments. 

If  the  admittance  were  known  as  a  function  of  pressure,  frequency  and 

^eatl^simnlifr^10?^  the,prediction  of  combustor  stability  would  be 
greatly  simplified.  There  is  no  sound  basis  for  such  predictions  from 

m°dels\  The  PurPose  of  the  present  work  is  to  determine  the 
mittance  or  related  functions  from  experimentally  observed  combustor 

the  Idm^tU  08  ^  ,T‘bwrner-  In  the  f°ll0win^  relation  between 

the admittance  and  the  burner  behavior  is  presented  as  a  basis  for  ex¬ 
plaining  the  data  reduction  procedure. 

T-burnprer^ng  JUSt  outside  the  (thin)  combustion  zone  in  the 

,  l  am^litudeB  velocity  and  pressure  oscillations  are  de¬ 

noted  by  (AU)m  and  AP,  and  the  phase  between  velocity  and  pressure  is 
denoted  by  6.  The  component  of  (AU)m  in  phase  with  pressure  is 

AU  *  (AU)m  cos  9  (13) 

SLrt;rts  per  unit  ttae  and  unit  —  °f 


E  -  h  (AP)  (AU)  (14) 

with  F  representing  the  same  quantity  as  in  Eq.  12.  This  term  can  also 
interpreted  as  the  work  done  by  a  pressure  AP  acting  on  a  unit  of 

urface  which  moves  a  distance  AU  each  unit  time  as  a  consequence  of  the 
pressure  action.  From  Eq.  14  it  follows  that  consequence  of  the 


9 


V 


E  =  h  (AP)2  (~)  -  -  h  (AP)2Y 
where  Y,  the  acoustic  admittance,  takes  the  form2 


Y  -  -  (|jf)  (16) 

In  terms  of  the  combustion-sound  field  interaction,  the  requirement  for 
amplification  is  that  E  has  to  be  positive  or  Y  negative.  For  practical 
purposes,  however,  the  minus  sign  in  Eq.  16  can  be  ignored. 

Propellant  Response  Function 


Equation  16  can  be  transformed  into  a  more  workable  expression  by 
applying  the  steady  state  conservation  equation 


Up  ■  rpg  ■  m 


where  U  is  the  mean  velocity  of  the  combustion  gases,  r  is  the  propel¬ 
lant  steady  state  burning  rate  at  mean  pressure,  ps  is  the  propellant 
density,  and  m  is  the  propellant  mass  burning  rate  per  unit  area. 
Assuming  that  Eq.  17  remains  valid  when  a  variation  in  U  equal  to  AU  is 
imposed,  It  follows  that 


AU  =  A(— )  -  (M  _  m  £§■) 
P  P  P* 


where  Am  and  Ap  are  the  corresponding  variations  in  mass  burning  rate 
and  gas  density.3  By  rearranging  and  using  Eq.  17,  Eq.  18  becomes 


AU  «=  U  (—  -  ^-) 

m  p 


Defining  the  fractional  variations  in  mass  burning  rate,  pressure, 
and  density  as 


Am  Ar 
m  r  "  V 


2  This  is  the  real  part  of  the  complex  acoustic  admittance  used  in 
many  references. 

3  Incremental  quantities  refer  to  component  in  phase  with  pressure. 


the  admittance  function  in  Eq.  17  can  then  be  written  as 


Y 


U  ,y  £\ 
p  'e  ~ 


(23) 


where  the  ratio  y/e  has  become  known  as  the  "pressure-coupled  propellant 
response  function",  and  represents  the  ratio  of  the  fractional  variation 
in  burning  rate  to  the  fractional  variation  in  pressure.  Essentially, 
this  is  described  by 


£  _  1*1 

e  r  AP 


(24) 


where  P/r  is  the  steady  state  and  Ar/AP  the  oscillatory  term  of  the 
response  function. 

Acoustic  studies  of  T-burner  cavities  (Ref.  6)  have  shown  that  the 
admittance  function  given  by  Eq.  23  requires  an  additional  term  in  order 
to  account  for  the  effect  of  mean  flow  in  the  burner,  which  is  caused 
by  the  combustion  gases  as  they  proceed  downstream  in  both  directions 
from  the  burner  ends.  The  effect  of  the  mean  flow  is  a  convective 
transport  of  acoustic  energy,  and  the  value  of  the  correcting  term  is 
given  by  -  (U/P)  (1/y);  therefore,  a  new  expression  for  an  equivalent 
"no-mean-flow  admittance"  is 


Y' 


U 

P 


(25) 


It  is  apparent  that  the  mean  flow  has  a  destabilizing  effect  (same  sign 
as  y/e).  It  should  be  noted,  however,  that  Eq.  25  is  no  longer  the 
admittance  function  of  a  surface,  as  is  the  case  with  Eq.  23,  the  re¬ 
lation  between  the  admittance  and  Y’  being 


Y' 


Y 


U  1 
P  y 


(26) 


t 

Sometimes  a  reduced  specific  admittance  function  is  defined  by  removing 
the  U/P  term,  i.e., 
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y'  =  Y'/(U/P)  =  y  -  ^ 


(27) 


At  high  frequencies,  when  isentropic  or  adiabatic  conditions  should  pre¬ 
vail  in  the  burner  cavity,  the  ratio  a/ £  exactly  matches  1/y,  and  the 
admittance  function  takes  the  form 


Y 


HF 


U  p 
P  e 


(28) 


Conversely,  when  conditions  become  isothermal  at  low  frequencies,  the 
ratio  o/e  acquires  value  unity,  and 


Y 


LF 


_  U  (H  _  111) 

P  ^e  y  } 


(29) 


The  next  step  is  to  determine  the  relationship  between  a  and  Y  and, 
through  it,  the  relationship  between  a  and  p/e.  On  one  hand,  the  acous¬ 
tic  energy  rate  in  terms  of  the  standing  wave  equation  is  given  by 

(AP)  2L 

E  =  2aE  =  2a  (30) 


On  the  other  hand,  considering  a  double-ended  T-burner  of  length  L, 
where  there  are  two  unit  areas  of  combustion  zone  per  unit  cross  section 
area  of  cavity,  the  acoustic  energy  rate,  according  to  Eq.  15,  will  be1* 


E  =  -  2  [h  (AP)2Y] 


(31) 


thus,  from  the  combination  of  Eq.  30  and  31,  it  follows  that 


Y 


aL 

2pc2 


(32) 


Introducing  the  relation  for  the  oscillation  frequency  of  the  first 
acoustic  longitudinal  mode  of  the  cavity 


f 


(33) 


4  (AP)S  in  Eq.  30  is  the  same  as  AP  in  Eq.  31  because  the  propel¬ 
lant  is  located  at  the  pressure  antinode. 
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Eq.  32  may  be  written 


Y 


ct 

4cpf 


(34) 


this  being  the  same  result  obtained  hv  -  A-tec 

Alternatively,  if  the  isentroDir  a™  different  method  in  Ref.  7. 

7  isentropic  acoustic  relationship  is  used,  namely. 


,2  = 

P 


(35) 


Eq.  34  becomes 


(36) 

2E!  o^P^L^rinffre^cy"0^”  SS  "J  S  ™rl- 

corrected  expression  to  account'for  Ihl  ^  admlttance  function  the 

T-burner,  Eq.  25,  it  ^  ?f5ect  °f  mean  flow  ln  the 

q  it  follows,  together  with  Eq.  17  and  34,  that 


4cps  f  f  <7  -  £  +  -) 
s  r  Gey 


(37) 


As  the  value  of  a  and  that  of  4coc  r/P  ana  r  ^ 

Independently,  it  has  become  customarv  tn  pv  u™™  °r  deta™ined 

of  p/e.  When  rearranged,  the  equation  fnr  V,Tess  the  results  in  terms 
the  form  q  for  w/e  at  high  frequencies  takes 


(38) 


factor  (o,/f)fy°whi2h  is  d^ripSe^rtto/r'in”  e<,uatlon  involves  one 
second  factor  <P/cpsr)  ,  which  ij  desfrlptlveS SS  a'"'10'’  and  * 
tion  of  the  propellant.  T-burnpr  no  a  stea^y  state  combus- 

“1  +  -2)  and  other  necessary  d"ta  Scent  ^  of  “  <as 

value  of  p/e  is  readily  calculated  The  n  5^°UrS6’  Ps)  from  whlch  the 
that  define  the  neutral  boundary  are  deter^inedT*8/01  01  ®qual  t0  Zero 
6  last  unstable  and  first  stable  tests  in  a  series?^^013^00  between 

<■/«.  Partly  becaus^the  *  average  TJ  T  ^  of 

thr-h  “  —  --  Serfthr^br"^She 
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combustion  zone  and  partly  because  of  the  assumption  that 


The  correction  for  velocity  of  sound  amounts  to  (c/c*)2,  where  c  repre¬ 
sents  the  average  speed  of  sound  in  the  cavity  and  c'  the  speed  of 
sound  at  the  combustion  zone.  Correction  for  the  heat  loss  effect  is 
discussed  in  Ref.  8,  but  was  not  used  in  the  present  work.  At  high  fre¬ 
quency,  when  the  thermal  losses  are  generally  small,  both  speeds  of 
sound  can  be  assumed  to  be  the  same  in  good  approximation. 

From  Eq.  37  it  can  be  seen  that  low  values  of  c,  r,  and  y/e  and 
high  values  of  P  are  conducive  to  stability  (making  due  allowance  for 
interdependence  among  these  variables) . 


14 


NWC  TP  4393 


T-BURNER  APPARATUS 


Burner  Assembly 

The  combustion  chamber  coikI nf  . 

jh, ir  °f  1 1/2,lnches  -4  - ^°“^s«^itt“1„fpipe 

propeuant  s^"°m^ 

IT*  TTTTIT'  XT'  S  “ll0D 

straight,  5-inch-long  nozzle  neck  hSJSa  anT^  5?  3  SUrge  tank  by  a 
and  an  outside  diameter  of  1  1/2  inches8  The  diameter  of  5/8  inch 

have  a  total  volume  of  about  6.5  ft^or’ll  200  t3nk  containers 

schematic  layout  of  the  burner  u Jr  ]  *1  ‘  F±gUre  2  shows  a 

work,  propellants  were  *-«„*- c  nstrumentation.  In  the  present 

cps  and  at  pressures  to  2  400  Ztl  in  *****  fr0m  5°°  t0  10’000 

steps.  The  frequen  v  of  8  convenle"t  pressure  and  frequency 

tudinal  mode ofdhe  cavL 3t  the  flrSt  W- 
burner;  the  mean  pressure  was  varieJTthrouvhddM  len8th  of  the 
the  burner  and  aurgaStaU„rk  ass'eS^  thr°Ugh  lnltlal  of 

standardized  as  shou^ln'l^blc0!  ^The'IT108’  15  burner  settings  were 
given  for  halfway  through  burn!™  .Tf  g3S  COlU"n  len*ths  b 
(3/4  inch),  the  pipe  section^  and  ?h  6  the  central  nozale  bl»ck 

(1/4  inch),  when  required  The  f 6  connectors  between  sections 

on  the  basis  of  an  average  speed  of  an  ^raquencies  %  "ore  calculated 
The  burner  volume,  excluding^he  neck  "is  denTi/t8  vlth  f reqi“!ncy ' 

amountpOfgPropaiianr  at  2,0000^  aboS?  TuoTmlTlTJT  T 

referl:hpereCsstu™e1e™Le-a8a)e200S'4onCh  "a8  ab  bb-e 

and  1,600  psig,  depending ^  upon’ the  TtailTT  ^  0r  <2)  2°°'  800’ 
the  propellant  in  the  ,7  “notable  pressure  range  exhibited  by 

frequencies  were  chosen  as  referenced’  Elevan  of  the  15  standardized 
In  addition,  selected  Lsts  wfra  n”d  fqaenCleS  f°r  each  P™«<r.  level, 
of  the  neutral  boundary  JncSne  ^  t0  delI'"it  the  upper  reaches 
d it ions.  includin8  the  upper  pressure  and  frequency  con- 
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TABLE  1.  Burner  Settings  and  Nominal  Frequencies 


Setting 

No. 

Section  Lengths 
in. 

L 

in. 

1/2  L 
in.-1 

fn 

cps 

c 

in/sec 

VB 

in3 

1 

3/4 

! 

3/4 

2.0 

0.250 

10,000 

40,000 

3.5 

2 

1.0 

1.0 

2.5 

0.200 

7,900 

39,500 

4.4 

3 

1.5 

1.5 

3.5 

0.143 

5,590 

39,000 

6.2 

4 

2.0 

2.0 

4.5 

0.111 

4,280 

38,500 

7.9 

5 

2.5 

2.5 

5.5 

0.0909 

3,450 

38,000 

9.7 

6 

3.0 

3.0 

6.5 

0.0769 

2,880 

37,500 

11.5 

7 

3.5 

3.5 

7.5 

0.0667 

2,460 

37,000 

13.2 

8 

4.0 

4.0 

8.5 

0.0588 

2,150 

36,500 

15.0 

9 

4.5 

4.5 

9.5 

0.0526 

1,890 

36,000 

16.8 

10 

6.0 

6.0 

12.5 

0.0400 

1,400 

35,000 

22.0 

11 

2.5 

2.5 

16.0 

0.0312 

1,080 

34,500 

28.2 

12 

3.5 

6.0 

6.0 

3.5 

20.0 

0.0250 

840 

33,500 

35.3 

13 

2.5 

10 

9.0 

3.5 

26.0 

0.0192 

620 

32,500 

45.9 

14 

3.0 

14 

11 

6.0 

35.0 

0.0143 

450 

31,500 

62.0 

15 

3.0 

20 

17 

6.0 

47.0 

0.0106 

320 

30,000 

83.0 

Instrumentation 

The  pressure  oscillations  in  the  burner  cavity  are  sensed  by  a  water- 
cooled  pressure  transducer  in  contact  with  the  back  face  of  one  of  the 
propellant  specimens.  The  transducer  is  located  1/6  diameter  off  center, 
and,  provided  that  the  thickness  of  the  propellant  disk  is  kept  to  about 
1/4  inch,  there  is  no  apparent  reduction  in  signal  strength  or  quality. 

The  signal  is  received  by  an  oscillator-demodulator  and  split  into  two 
channels.  One  channel  passes  through  a  low-pass  DC  filter  and  is  dis¬ 
played  on  an  oscilloscope,  and  the  other  channel  is  taken  by  a  high-pass 
AC  filter  and  again  split  into  two  channels.  One  channel  is  then  dis¬ 
played  as  an  AC  trace  on  the  same  oscilloscope  and  the  other  is  directed 
to  a  low-pass  AC  filter.  The  resulting  signal  is  finally  amplified  and 
recorded  on  photographic  chart  paper  by  means  of  a  galvanometer  oscillo¬ 
graph.  The  oscilloscope  displays  are  then  recorded  on  Polaroid  film. 

The  high-  and  low-pass  AC  filters  ensure  that  the  signal  registered  by 
the  galvanometer  oscillograph  contains  only  the  first  longitudinal  mode 
of  the  acoustic  pressure  oscillations  in  the  burner  cavity. 

The  oscilloscope  record  provides  an  immediate  visualization  of  the 
test  and  simultaneously  presents  the  mean  burner  pressure  and  the  oscil¬ 
lation  pressure.  Since  the  higher  frequencies,  if  present,  are  unfiltered 
at  this  point,  the  record  also  indicates  the  content  of  higher  oscillation 
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modes.  The  mean  burner  pressure  trace  is  used  to  calculate  the  burning 
rate  of  the  propellant  at  that  pressure.  The  oscillograph  record  is  the 
main  source  of  data  acquisition  and  measures  the  actual  frequency,  04  and 
(*2»  and  the  maximum  peak-to-peak  value  of  the  acoustic  pressure  amplitude, 
(AP)m.  The  record  provides  an  indication  of  the  quality  of  the  test, 
showing  the  effect  of  the  igniter  pulse,  burning  symmetry,  and  distor¬ 
tions  at  burnout,  from  which  it  can  be  decided  whether  or  not  to  repeat 
the  test. 

At  least  three  calibrations  are  needed  to  ascertain  instrumentation 
performance:  (1)  static  pressure  calibration  of  the  pressure  transducer 

by  means  of  a  deadweight  testing  machine  at  various  attenuations  in  the 
oscillator-demodulator  unit,  (2)  calibration  of  the  frequency  response  of 
the  pressure  transducer  and  AC  filter  units  at  constant  voltage,  and  (3) 
calibration  of  the  amplitude  response  of  the  pressure  transducer  and  AC 
filter  units  at  constant  frequency.  Calibration  (1)  is  carried  out 
periodically  to  check  deterioration  or  malfunction  of  the  transducer 
probe,  and  calibration  (2)  is  made  whenever  the  AC  filter  units  are  sus¬ 
pected  to  be  faulty.  The  pressure  transducer  response  is  almost  inde¬ 
pendent  of  frequency  up  to  about  8,000  cps ,  but  falls  sharply  at  in¬ 
creasingly  higher  frequencies. 

Test  Procedure 

The  test  procedure  comprises  the  following  steps: 

1.  Machining  the  propellant  specimens  to  disk  shape  with  a  nominal 
thickness  of  1/4  inch  and  a  diameter  of  1  1/2  inches 

2.  Fitting  and  bonding  the  disks  to  the  terminal  burner  sections 
with  epoxy  resin 

3.  Curing  the  epoxy  bond 

4.  Fixing  the  igniter  (spreading  pyrotechnic  material  over  the 
surface  of  the  propellant  disk  and  teardrop  primer  of  the  same 
substance,  which  supports  the  ignition  wires 

5.  Assembling  the  burner  sections  and  end  caps 

6.  Soldering  the  ignition  wires  inside  the  burner  cavity  to  the 
terminal  posts  at  the  nozzle  neck 

7.  Mounting  the  burner  unit  to  the  surge  tank  and  connecting  the 
transducer  probe 

8.  Pressurizing  the  burner  atmosphere  with  nitrogen  to  the  desired 
pressure  level 

9.  Checking  on  instrumentation  settings  and  circuit  continuity 
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10.  Activating  the  sequencer,  igniting,  and  recording  the  test 
data. 

Data  Acquisition 

The  data  acquisition  procedure  is  aimed  at  determining  all  the 
necessary  terms  involved  in  Eq.  38.  The  steady  state  terms  are  readily 
known  or  calculated.  P  is  the  initial  pressurization  level  imposed  on 
the  burner,  r  is  calculated  from  the  oscilloscope  DC  pressure  trace  or 
from  strand  burning  rate  measurements,  and  the  value  of  ps  is  easily 
obtained. 

The  oscillatory  terms  are  derived  from  the  oscillograph  record. 

The  actual  frequency  of  the  pressure  oscillations  is  read  after  it  has 
attained  a  constant  value  during  the  growth  of  the  oscillations,  and  c 
is  determined  as  the  average  slope  in  the  plot  of  f  versus  1/2  L,  L 
being  the  length  of  the  burner  gas  column  halfway  through  burning  (see 
Table  1)  .  Most  of  the  data-reduction  work  goes  into  determining  the 
growth  and  decay  constants  of  the  pressure  oscillations.  A  standard 
numerical  procedure  renders  the  average  value  of  aj ;  the  value  of  f 
corresponds  to  this  stage  in  the  test,  and  in  most  cases  it  remains  sub¬ 
stantially  the  same  until  the  moment  of  burnout.  The  calculation  of 
is  not  so  straightforward  and  implies  a  certain  degree  of  arbitrari¬ 
ness.  The  conditions  at  burnout  are  such,  mainly  through  thermal  losses, 
that  the  decay  constant  of  the  pressure  oscillations  usually  undergoes  a 
noticeable  change  from  beginning  to  end  of  the  decay  period;  the  tendency 
is  toward  smaller  values  (although  sometimes  the  opposite  is  true  at  high 
frequencies)  and  toward  lower  corresponding  frequencies  as  the  decay  pro¬ 
ceeds  in  time.  Consequently,  each  decay  event  is  represented  by  two 
constant  and  frequency  combinations,  one  denoted  by  (012,  f^) *  reflecting 
the  beginning,  and  the  other  by  (02,  f2) »  reflecting  the  end  of  the  decay 
period.  The  maximum  peak-to-peak  value  of  the  amplitude  of  the  acoustic 
pressure  oscillations  (AP)m,  is  calculated  from  the  records  with  the  help 
of  the  calibrations  mentioned  above. 

The  determination  of  the  final  value  of  a  comprises  a  sequence  of 
plotting  steps  intended  to  smooth  the  raw  test  data.  The  following  steps 
have  been  adopted:  (1)  plots  of  aj  versus  f  and  02  versus  f2,  whereby 
each  test  is  represented  by  a  line  extending  from  (02,  £2)  to  (02,  f^') 
in  the  latter  plot,  (2)  graphical  averaging  of  the  data  by  a  best  fit 
curve,  (3)  tabulation  of  the  04  and  c*2  values  for  the  reference  fre¬ 
quencies  500,  1,000,  2,000  cps,  etc.,  (4)  replots  of  these  data  as  a^/f 
versus  f  and  c^/f  versus  f  and  further  graphical  averaging  if  necessary. 
The  values  of  a  are  given  by  the  final  04  and  02  values,  and  the  pro¬ 
cessed  data  are  collected  in  tabular  form  as  presented  under  "Experimen¬ 
tal  Results."  The  value  of  the  limiting  frequencies  for  04  =  0  (neutral 
boundary  frequencies)  are  also  determined  graphically  from  the  a*  versus 
f  plots  for  each  pressure. 
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Attempts  were  made  to  obtain  directly  the  logarithmic  recordings 
of  the  growth  and  decay  of  the  pressure  oscillations  by  incorporating  a 
Tektronix  (Type  0)  logarithmic  converter  amplifier  into  the  instrumen¬ 
tation.  Had  this  approach  been  successful,  the  very  time-consuming 
numerical  procedure  for  determining  aj  and  012  would  have  been  greatly 
reduced . 


Accuracy  of  Results 

There  is  a  multitude  of  possible  sources  of  error  that  can  adverse¬ 
ly  affect  the  quality  of  the  results.  They  fall  into  three  major  cate¬ 
gories:  (1)  the  nature  of  the  technique  under  which  the  T-burner  nystem 

operates,  (2)  the  measuring  equipment  which  senses  and  records  the  pri¬ 
mary  data,  and  (3)  the  method  of  data  acquisition. 

The  error  contributions  by  the  technique  concern 

1.  Propellant  manufacture,  whereby  batch-to-batch  variations 

in  composition  or  quality  tend  to  impair  test  repeatability, 
especially  as  marginally  unstable  conditions  are  approached 
near  the  neutral  boundary 

2.  Burner  atmosphere,  initially  composed  of  almost  pure  nitrogen, 
but  becoming  richer  in  combustion  products  as  the  propellant 
burns,  with  the  result  that  the  growth  and  decay  of  the  pres¬ 
sure  oscillations  must  occur  in  an  environment  of  changing 
temperature  and  composition,  which  weakens  the  hypothesis  of 
Constance  of  attenuation,  especially  for  propellants  that 
give  off  a  considerable  amount  of  particulate  matter 

3.  Igniter  pulse,  produced  by  the  rapid  burning  of  the  pyrotechnic 
material,  which  can  distort  the  growth  profile  of  the  pressure 
oscillations,  and  also  the  presence  of  the  ignition  wires  and 
their  combustion  products  inside  the  burner  cavity,  since 
their  effect  as  acoustic  attenuators  cannot  be  ruled  out 

4.  Nonsimultaneous  or  nonuniform  burning  of  the  propellant  speci¬ 
mens  caused  by  different  ignition  delays  or  failure  of  the 
propellant-to-metal-wall  bonding 

5.  Variations  in  operating  conditions,  such  as  initial  propel¬ 
lant  temperature  or  mean  burner  pressure,  when  the  initial 
level  of  pressurization  cannot  be  maintained  in  the  surge  tank. 

In  general,  these  sources  of  error  become  more  critical  as  low  fre¬ 
quencies  are  approached.  This  is  particularly  true  of  (1)  the  burner 
atmosphere  (because  of  the  increasing  larger  burner  volume  as  the  fre¬ 
quency  is  reduced,  which  gives  rise  to  attempts  to  arrest  the  growth  of 
the  pressure  oscillations  (hold-off  grains,  movable  acoustic  attenuators) 
until  sometime  later  during  burning)  and  (2)  the  igniter  pulse  (because 
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it  tends  to  persist  longer  at  low  frequencies  because  of  small  cavity 
attenuation  and  because  of  the  long  length  of  ignition  wire  present  in 
the  cavity  at  low  frequencies,  much  of  which  remains  unconsumed).  These 
combined  effects  have  limited  extention  of  routine  testing  below  about 
500  cps,  which  is  considered  to  be  the  lowest  practical  limit  of  the 
present  T-burner.  Uncertainties  caused  by  propellant  composition  or 
quality  are  usually  small;  as  for  the  unsymmetry  in  the  burning  of  the 
propellant  specimens,  which  can  happen  at  any  frequency,  record  dis¬ 
tortion  is  normally  great  enough  for  it  to  pass  undetected. 

Instrumentation  malfunctions  can  cause  errors  in  oscilloscope  and 
galvanometer  oscillograph  records.  The  three  components  most  liable  to 
malfunction  are  the  pressure  transducer  probe,  the  oscillator-demodulator, 
and  the  AC  filters.  The  pressure  transducer  tan  be  recalibrated  periodi¬ 
cally  or  whenever  visual  inspection  shows  a  sudden  deterioration  of  the 
sensing  head.  Marginal  malfunctions  of  the  demodulator-oscillator  and 
AC  filters  can  go  undetected  for  several  tests  until  the  warning  symp¬ 
toms  become  severe.  The  appearance  of  minor  spurious  features  in  the 
records  is  not  uncommon  and  usually  can  be  disregarded  unless  they  be¬ 
come  repetitive. 

The  data  reduction  procedure,  though  simple  in  itself,  depends  on 
the  operator's  judgment,  particularly  for  poorly  defined  records.  In 
general,  the  array  of  uncertainties  is  such  that  no  estimation  of  abso¬ 
lute  errors  is  possible;  the  only  way  to  assure  accurate  results  is  by 
checking  repeatability  and  by  adopting  a  fixed  procedure  for  testing, 
data  acquisition,  and  data  reduction.  Scattering  of  test  data  points 
indicates  that  the  final  and  response  function  values  have  an  uncertainty 
no  better  than  +10%  when  04  and  02  are  both  above  10  sec-1  and  provided 
04  is  not  much  smaller  than  02*  This  is  generally  the  case  for  condi¬ 
tions  of  medium  to  severe  instability  in  the  high-frequency  region. 

At  low  frequencies,  when  04  and  02  might  fall  well  below  10  sec-* ,  or 
near  the  high-frequency  side  of  the  neutral  boundary,  when  o  approaches 
a2*  the  uncertainty  increases  to  +15%  or  higher. 
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PROPELLANTS 


Selection  of  Propellants 

Most  of  the  propellants  tested  belong  to  an  "A"  series  of  research 
composite  propellants  formulated  to  encompass  a  wide  number  of  composi¬ 
tional  variations;  others  are  commercial  propellants,  including  four 
typical  double-base  propellants  developed  years  ago.  Among  the  38  pro¬ 
pellants  studied,  extended  testing  was  confined  to  nine  composite  formu¬ 
lations  (A-13,  A- 14 ,  A-15,  A-16,  A-17,  A-18,  A-19,  A-81,  and  T-35)  and 
four  double-base  formulations  (N-4,  N-5,  JPN,  and  X-14) .  This  selection 
was  made  in  order  to  characterize  some  of  the  major  compositional  and 
burning  rate  effects  on  instability,  such  as  (1)  differences  of  behavior 
between  composite  and  double-base  propellants,  (2)  influence  of  the 
linear,  mesa,  and  plateau  burning  rate  modalities  for  the  double-base 
propellants,  (3)  effect  of  ammonium  perchlorate  particle  size  and  com¬ 
bustion  catalysts,  and  (4)  effect  of  binder,  oxidizer  types,  and  aluminum 
for  the  composite  propellants.  Of  the  remaining  propellants,  only 
partial  evaluation  is  available  for  16  (A-l,  A-2,  A-5,  A-7,  A-10,  A-ll, 
A-12,  A-20,  A-22,  A-23,  A-24,  A-25,  A-26,  A-27,  A-28,  and  A-44) ,  and 
five  propellants  (A-52,  A-83,  A-91,  A-138,  and  A-139)  did  not  exhibit 
instability  within  the  frequency  covered  by  the  present  T-burner. 

Propellant  Procurement 

The  A-series  propellants  are  manufactured  by  the  Development  Engi¬ 
neering  Branch,  Propulsion  Development  Department,  at  NWC,  using  routine 
mixing  and  curing  procedures.  Of  the  three  binders  used,  PBAN,  PU,  and 
PS,  good  quality  is  observed  for  the  propellants  containing  PBAN,  only 
marginal  quality  for  those  containing  PU,  and  attempts  to  produce  PS 
formulations  resulted  in  propellants  of  unacceptable  quality.  Some 
quality  alteration  has  been  noticed  for  the  PU  binder  propellants  in 
that  aging  at  room  temperature  seems  to  bring  about  a  tendency  toward 
stability. 

Each  batch  yielded  approximately  7  pounds  of  propellant  or  about 
4  pounds  of  machined  propellant  specimens,  providing  material  for  80 
tests,  a  number  ample  enough  to  evaluate  some  of  the  less  unstable 
propellants. 
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Propellant  Composition  and  Burning  Rate 

Composite  compositions  are  shown  in  Table  2,  and  double-base  com¬ 
positions  are  shown  in  Table  3.  The  burning  rates  are  shown  graphically 
in  Fig.  3  and  4  for  the  composite  and  in  Fig.  5  for  the  double-base  pro¬ 
pellants.  The  composite  propellants  and  the  N-4  and  JPN  double-base 
propellants  exhibited  essentially  linear  burning  rates  in  the  usual 
logarithmic  plot  of  r  versus  P;  the  N-5  propellant  exhibited  mesa,  and 
the  X-14  plateau  burning  characteristics.  It  should  be  noted  that  the 
PBAN  binder  ingredient  in  the  A  propellants  has  been  mistakenly  reported 
in  the  past  as  being  PBAA.  The  A-28  propellant  was  not  formed  by  cast¬ 
ing  but  through  pressing  the  mixture  at  40,000  lb/in2. 

The  85y  AP  was  a  commercial  grade  designated  Lot  90-3  and  90-4  by 
the  American  Potash  and  Chemical  Corporation,  and  the  15y  AP  was  stan¬ 
dard  American  Potash  and  Chemical  Corporation  15p  grade.  Ammonium  per¬ 
chlorate  designated  by  BM  in  the  tables  and  figures  was  a  50-50  mixture 
of  micropulverized  (approximately  9p)  and  as-received  (approximately 
150y)  material. 


I 


A 


mm m 


TABLE  2.  COMPOSITE  PROPELLANTS 
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Binder:  /„  wt.  includes  curing  agents;  rBAN  =  polybutylacrylonitrile  acrylic  acid,  PS  =  polv- 
sulphide,  PU  =  polyurethane.  Oxidizer:  AP  =  ammonium  perchlorate,  KP  -  potassium  perchlorate, 

BM  =  bimodal  distribution  containing  micropulverized  and  “'as  received"  material  in  equal  parts. 
Catalysts:  cc  =  copper  chromite,  LF  =  lithium  fluoride,  CB  =  carbon  black. 


TABLE  3(a).  Double-Base  Propellants 
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Nitrocellulose 
Nitroglycerin 
Diethl  phthalate 
Ethyl  centralite 
Carbon  black 
Potassium  sulfate 
Candelilla  wax 
2-Nit rodiphenylanime 
Lead  stearate 
Lead  salicylate 
Lead  2-ethylhexoate 


JPN 

N-4 

N-5 

51.40 

51.0 

50.0 

42.90 

34.3 

34.9 

3.23 

10.6 

10.5 

1.00 

— 

0.20 

0.1 

1.25 

1.5 

_ 

0.02 

— 

0.2 

—  — 

2.0 

2.0 

“  “ 

0.5 

—  — 

— 

1.2 

1.2 

100.0% 

100.0% 

Ioo.o% 

TABLE  3(b).  Burning  Characteristics 


JPN 

N-4 

1 - 

N-5 

P  »  lb/in3 . 

.059 

.056 

.056 

r»  in/sec  at  P,  psia 

215 . 

.22 

415 . 

.17 

.16 

815 . 

— — 

1615 . 

.55 

.34 

.44 

.89 

.48 

.41 

Burn  rate  equation 

r  =  CPn,  in/sec  .  .  .  . 

5.4X10'V69 

11.0X10~V51 

Mesa  I 

EXPERIMENTAL  RESULTS 


General  Form  of  Behavior 


When  a  T-burner  is  fired  under  unstable  conditions,  oscillations 
develop  shortly  after  ignition,  grow  exponentially,  level  off  at  some 
steady  amplitude,  and  decay  more  or  less  exponentially  after  burnout. 
Most  unmetallized  propellants  produce  unstable  behavior  under  some  con¬ 
ditions  of  pressure  and  frequency,  and  determinations  can  thus  be  made 
of  the  combustion  a  or  response  function  (real  part)  y/e.  Tests  are 
fired  at  a  specific  mean  pressure  with  burners  of  several  lengths,  and 
the  series  is  repeated  for  other  mean  pressures  as  needed.  From  such 
tests,  a  graph  of  response  function  versus  frequency  is  developed,  with 
curves  for  various  mean  pressures  (Fig.  6).  In  an  extensive  series  of 
tests,  the  range  of  burner  lengths  (frequencies)  tested  will  usually 
include  length  extremes  for  which  the  burner  is  stable.  These  sets  of 
tests  delineate  not  only  the  response  characteristics  of  the  propellant, 
but  also  the  stability  limits  of  the  propellant-combustor  combination. 
These  combined  properties  may  be  summarized  in  a  diagram  such  as  Fig.  7, 
which  is  a  contour  map  of  the  response  function  in  the  pressure- 
frequency  plane.  The  outer  boundaries  are  the  stability  limit  for  the 
type  of  burner  used,  and  iso-response  contours  obtained  by  cross  plot¬ 
ting  in  the  y/t  versus  f  diagram  represent  the  perturbation  combustion 
characteristics  of  the  propellant.  The  contour  map  shows  in  this  case 
a  maximum  response  function  at  point  P*f*.  The  two  curves  through  that 
point  are  the  locus  of  frequencies  for  maximum  y/e  as  a  function  of 
mean  pressure,  and  locus  of  pressures  for  maximum  y/c  as  a  function  of 
frequency  (e.g.,  the  ff  curve  in  Fig.  7  is  a  plot  of  the  maximum  y/e 
points  of  the  curves  in  Fig.  6).  Ideally,  one  would  like  to  have  a 
three-dimensional  plot  of  these  functions  for  each  propellant  of 
interest,  but  cost  considerations  will  usually  limit  propellant  charac¬ 
terization  to  determination  of  the  high  y/e  portions  of  curves  like 
those  in  Fig.  6 


Response  Function 


The  response  function  values  have  been  calculated  from  Eq.  39,  i.e.. 


I*b)( 


cti  +  c*2 


and  a  being  equal  to  04  +  <*2,  an  abbreviation  for  each  pressure  level 
can  be  written  as 
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where  Kr  is  the  shortened  form  of  P/4cpsr  and  is  dimensionless. 

The  response  function  results  are  separated  into  two  groups  accord¬ 
ing  to  the  pressure  coverage. 

Group  I  comprises  the  results  obtained  testing  propellants  under 
limited  pressure  coverage,  mainly  at  200  psig  only;  these  results  are  to 
a  great  extent  recalculations  from  data  already  published  in  Ref.  9  and 
10,  together  with  additional  data  obtained  to  check  previous  results  or 
to  complete  missing  information.  These  results,  which  contain  only 
response  function  determinations  (since  no  other  information  can  be  re¬ 
liably  derived  from  scant  pressure  coverages)  are  presented  collectively 
in  Appendix  A  along  with  response  function  diagrams.  The  results  for 
the  A-44  propellant  are  suspect  because  of  the  poor  quality  of  this 
propellant. 

Group  II  includes  the  results  obtained  testing  propellants  under 
extensive  pressure  coverage;  previous  data  on  some  of  these  propellants, 
Ref.  9  and  11,  were  entirely  retested.  These  results  include  determin¬ 
ation  of  combustor  stability  limits  as  well  as  response  function.  The 
neutral  boundary  limits,  the  pressure  and  frequency  peak  lines,  and  the 
vertex  point  are  illustrated  in  the  P-f  diagrams,  also  in  Appendix  B, 
along  with  response  function  diagrams  for  the  three  corresponding  refer¬ 
ence  pressures.  The  conditions  yielding  maximum  response  function  were 
calculated  by  a  method  described  in  Appendix  C,  and  are  shown  in  Table  4 
for  those  propellants  that  were  tested  extensively. 


TABLE  4.  Vertex  Points  (Maximum  Response  Function) 


Propellant 
Serial  No. 


Propellant 
Serial  No. 
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Stability  Limits 

For  those  propellants  that  were  most  extensively  tested,  a  rough 
estimate  of  the  stability  limit  of  the  combustor  can  be  determined. 

This  is  illustrated  in  Fig.  8  for  a  hypothetical  example,  where  the 
value  of  04  is  plotted  as  a  function  of  frequency  for  a  particular  mean 
pressure.  When  the  curve  is  extrapolated  to  04  "0,  the  stability  limit 
is  given  by  the  corresponding  values  of  pressure  and  frequency. 

In  principle,  a  sufficient  number  of  oj  =  0  points  (p,  f)a=o  can  be 
obtained  in  this  way  to  plot  a  stability  limit  profile  like  that  in 
Fig.  7.  Unfortunately,  the  scatter  of  measurements  becomes  very  bad  as 
the  stability  limits  are  approached,  so  that  determination  of  the  limits 
is  accomplished  with  only  marginal  accuracy. 

In  order  to  achieve  a  more  systematic  correlation  of  the  rather 
limited  data  on  stability  limits,  a  computational  scheme  was  developed 
for  constructing  stability  boundaries  from  the  limit  points.  This 
scheme,  discussed  in  Appendix  C,  fits  the  limit  points  with  segments  of 
two  curves,  one  bounding  the  high  frequency  side  of  the  unstable  region, 
the  other  bounding  the  low  frequency  side.  These  segments  can  each  be 
fitted  when  three  limit  points  are  available,  and  can  then  be  described 
in  terms  of  limiting  high  and  low  values  of  pressure  and  frequency  for 
the  unstable  region.  Table  5  lists  these  values,  estimated  for  those 
propellants  having  sufficient  data  for  computations.  The  meaning  of 
the  notation  is  clarified  in  Fig.  7.  Table  4  listed  the  vertex  points 
(P*,f*)  for  those  propellants  with  extensive  test  data. 

Acoustic  Pressure  Results 

The  maximum  acoustic  pressure  levels,  (AP)m,  have  been  calculated 
for  most  of  the  propellants  of  the  extended  pressure  series.  These 
values  are  easily  derived  from  the  test  records  and  corresponding  cali¬ 
brations.  No  assessment  can  be  made  of  their  accuracy,  which,  on  the 
average,  should  not  be  regarded  as  very  high.  It  has  sometimes  been 
the  practice,  mainly  from  an  engineering  viewpoint,  to  describe  the 
degree  of  combustion  instability  in  terms  of  the  pressure  ratio  (AP)m/P, 
although  this  representation  has  not  yet  been  theoretically  substantiated. 
The  (AP)m  pressure  levels  and  (AP)m/P  pressure  ratios  for  the  propellants 
are  listed  in  Appendix  D.  Comparison  with  other  instability  data  seems 
to  suggest  a  certain  resiemblance  of  trends  between  the  04 /f  and  (AP)m/P 
ratios.  Also  relevant  !:o  this  large  amplitude  aspect  of  combustion 
instability  (but  not  explored  here)  is  the  change  on  propellant  burning 
rate  effected  by  the  amplitude  and  frequency  of  the  acoustic  pressure 
oscillations  sustained  in  the  burner  cavity,  as  the  experimental  ob¬ 
servations  in  Ref.  12  and  13  have  clearly  shown. 
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FIG.  8.  Plot  of  Versus  f  and  Extrapolation 
to  Determine  Stability  Limits. 


Effect  ol  Additives  on  Attenuation 


The  presence  of  particulate  material  in  the  reaction  products  is 
known  to  produce  a  major  effect  on  attenuation  of  oscillations.  Figure 
9  shows  the  trends  of  f°r  several  of  the  propellants  tested  on  this 
program.  A  comparison  of  the  results  for  a  propellant  with  "clean" 
reaction  products  such  as  A-14  with  the  other  propellants  shows  a  con¬ 
spicuous  effect  even  for  the  low  concentrations  of  condensed  products 
resulting  from  2%  copper  chromite  or  1%  aluminum. 

Detailed  Study  of  A-13  Propellant 

Considering  that  the  most  adequate  description  of  behavior  is 
given  by  the  P-f  diagrams,  showing  not  only  the  structural  lines  but 
also  the  isoresponse  function  contours  at  discrete  intervals,  further 
testing  with  the  A-13  propellant  was  undertaken  to  determine  some  of 
these  contours.  The  A-13  propellant  was  selected  because  (1)  it  has  a 
simple  formulation  and  is  easy  to  manufacture,  giving  batches  of  con¬ 
sistent  quality,  and  it  does  not  appreciably  deteriorate  with  time;  and 
(2)  its  unstable  behavior  is  confined  to  a  pressure  and  frequency  range 
that  is  attainable  in  the  present  T-burner  apparatus. 

The  testing  was  extended  in  two  directions:  (1)  toward  low  pres¬ 
sures  to  resolve  the  steep  gradient  of  the  isoresponse  function  contours 
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as  low  pressure  and  frequencies  are  approached,  and  (2)  at  constant- 
frequency  conditions  to  further  ascertain  contour  location.  A-13  growth 
and  decay  constants  are  shown  as  functions  of  frequency  in  Appendix  E. 
The  figures  contain  individual  test  data  points,  showing  the  scatter 
that  is  encountered.  In  cases  where  the  decay  assumes  two  values,  the 
two  points  are  connected  with  a  thin  straight  line.  The  assumed  vari¬ 
ations  of  growth  and  decay  constants  are  shown  as  heavier  curved  lines. 
The  cross-plotting  of  data  into  p-f  diagrams  has  produced  some  of  the 
isoresponse  function  contours  as  well  as  the  graphical  representation  of 
other  parameters  of  interest,  including  a,  a*,  and  «2  and  their  corre¬ 
sponding  frequency  ratios.  The  isoresponse  function  contours  and  the 
remaining  diagrams  are  given  in  Appendix  E.  The  specimen  results 
gathered  from  the  A-13  propellant  illustrate  additional  trends  that 
should  be  incorporated  into  the  general  form  of  behavior. 

High-Speed  Photography  of  Burning  Surfaces 

High-speed  cinematography  of  propellant  samples  burning  under 
steady  state  conditions  has  afforded  observation  of  surface  and  gas 
phase  events  that  have  some  qualitative  relevance  on  instability  be¬ 
havior.  The  observations  are  summarized  for  the  following  propellants. 

1*  N-4.  At  200  psig,  black  melt  and  moving  dark  particles  (spher¬ 

ical,  100  to  200  microns)  on  surface,  no  flame  visible;  at  400  psig, 
larger  number  of  dark  particles  smaller  in  size,  no  flame  visible;  at 
800  psig,  same  trend  for  particles,  appearance  of  flame  about  1,000 
microns  from  surface,  well  defined  dark  zone. 

2.  N-5.  At  200  psig,  irregular  particles  somewhat  spectacular 
(about  200  microns)  accumulate  on  surface  prior  to  shedding,  no  flame 
visible;  at  400  psig,  almost  same  appearance  as  before;  at  800  psig, 
larger  number  of  dark,  shiny  particles,  spherical  in  shape  as  they  move 
through  the  flame,  which  appear  about  1,500  microns  above  surface . 

3.  X-14.  At  400  psig,  multitude  of  small  shiny  blobs  on  surface, 
no  flame  visible;  at  800  psig,  visible  flame,  the  particles  appear 
spherical  and  incandescent  as  they  move  through  the  flame. 

Information  on  composite  propellants  is  summarized  below: 

1*  A-13.  At  200  psig,  dark  blobs  of  molten  binder  (150  microns) 

on  surface,  orange  flame  about  1,800  microns  high;  at  400  psig,  blobs 
of  molten  binder  tend  to  coalesce,  more  intense  and  uniform  flame  about 
1,500  microns  high;  at  800  psig,  blobs  of  molten  binder  increase  in 
number  and  size,  more  vigorous  and  higher  flame  than  before. 

2*  A- 14.  At  200  psig,  no  binder  melting,  weak  flame;  at  400  psig, 
no  melting,  stronger  flame;  at  800  psig,  no  melting,  stronger  flame, 
particles  eject  from  surface,  some  ammonium  perchlorate  crystals  appear 
to  protrude  as  surface  regresses. 
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3.  A-15.  Same  behavior  as  with  A-13,  but  more  intense  melting  of 

binder,  occasional  green  flares  appear  to  originate  from  the  burning 
binder  blobs. 


4.  A-16.  Same  behavior  as  with  A-14  but  more  vigorous  flame  and 
absence  of  ammonium  perchlorate  protrusion. 

5.  A-17 .  At  200  psig,  less  melting  of  binder  and  in  smaller  blobs 
than  with  A-13,  deposition  of  soot  on  surface,  weak  flame;  at  400  psig, 
considerable  increase  of  molten  binder  in  blobs  of  larger  size,  some  in 
filament  form,  stronger  flame;  at  800  psig,  same  trend  as  before,  but 
distinct  flares  appear  to  originate  from  burning  binder  blobs. 

6.  A- 18.  At  200  psig,  similar  to  A-14;  did  not  sustain  burning 
at  higher  pressures. 

7.  A-19.  At  200  psig,  no  binder  melting;  at  400  psig,  incipient 
binder  melt;  at  800  psig,  enhanced  melting  of  binder  and  intense  flame. 

8.  T-35 .  At  200  psig,  slight  binder  melt,  burning  metallic  parti¬ 
cles  collect  on  surface  and  later  move  through  flame  (resemblance  to 
aluminum  combustion);  at  400  and  800  psig,  same  behavior  but  with  fewer 
metal  particles  and  increasingly  stronger  flame. 

9.  A-52,  A-83,  A-91,  A-139.  Since  these  are  aluminized  propellants, 
the  most  prominently  observed  feature  is  the  combustion  of  the  metal  in¬ 
gredient;  aluminum  combustion  appears  to  enhance  melting  of  the  binder. 

10.  A-138.  Surface  almost  covered  by  a  bubbling  mass  of  molten 

potassium  perchlorate  and  very  strong  flame. 
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In  addition  to  the  stability  characteristics  of  the  propellants, 
there  is  reason  for  interest  in  the  stability  limits  of  the  T-burner. 

In  the  present  work,  considerable  emphasis  was  given  to  testing  near 
the  stability  limits,  and  on  correlation  of  results  by  analysis.  Re¬ 
sults  of  these  studies  should  be  of  value  in  future  development  of 
T-burner  testing  and  theory,  but  are  not  discussed  extensively  here. 

Effect  of  Binder  on  Response  Function 

The  test  series  included  a  comparison  of  three  propellants  that 
differed  only  in  the  type  of  binder.  The  response  functions  for  these 
propellants  (A12,  A22,  and  A25)  are  collected  in  Fig.  10.  The  response 
function  was  roughly  the  same  for  all  three  binders  at  the  low  end  of 
the  frequency  range  tested,  but  differed  significantly  at  high  frequency, 
with  polyurethane  being  most  unstable  and  PBAN  being  most  stable.  Mak¬ 
ing  allowances  for  the  different  values  of  r  for  these  propellants 
(PBAN,  r  =  0.21  in/sec;  PS,  r  =  0.29  in/sec;  PU,  r  *  0.19  in/sec),  sta¬ 
bilities  are  roughly  equal  in  the  high  frequency  range,  with  PS  propel¬ 
lant  being  the  most  unstable  at  "low"  frequency.  This  generalization 
should  not,  of  course,  be  extended  to  all  pressures  or  other  propellant 
formulations  without  considering  trends  such  as  those  discussed  below. 


FIG.  10.  Effect  of  Binder  on  Response  Function  (80%  AP,  50-50 
Mix  of  Micropulverized  and  As-Received;  No  Ballistic  Modifier). 
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Effect  of  Equivalence  Ratio 

AP  JriLT  made  0n  a  family  0f  Pr°Pellants  with  PBAN  binder  and  80w 
;*Wlth  °Xidi*er  contents  Of  76,  85  and  90%  AP  (A13,  A27,  A28) .  The 

valuesSof  rCt!°nS  316  Sh0Wn  in  Fi8‘  11  f°r  tests  at  215  Psia,  along  with 
values  of  burning  rate.  From  the  figure  it  is  evident  that  the  values 

the  response  functions  differ  conspicuously,  although  they  usually 

£rr 


FIG.  11. 
Function 
AP). 


Effect  of  Equivalence  Ratio  on  Response 
(PBAN  Binder,  No  Ballistic  Modifier,  80y 


This  variable  was  examined  using  a  PBAN  binder,  an  oxidizer  content 
of  76%,  and  using  either  lithium  fluoride  (LiF) ,  copper  chromite  (cc) , 
or  no  ballistic  modifier.  Results  are  summarized  in  Fig.  12(a)-(d), 
which  differ  primarily  according  to  ballistic  modifier.  The  formula¬ 
tions  with  no  modifier  (Fig.  12(a)  and  (b))  were  tested  at  215  psia  and 
815  psia,  with  15y,  bimodal  (50-50  mix  of  15y  and  80y) ,  and  80y  particle 
sizes.  The  formulations  with  ballistic  modifiers  were  tested  with  only 
15y  and  80y  unimodal  oxidizer  (Fig.  12(c)  and  (d)).  The  values  of  y/e 
were  better  correlated  by  f/r2  than  by  f,  and  were  thus  plotted  in  these 
figures  as  y/e  versus  f/r2. 

In  all  cases  where  the  oxidizer  particle  size  was  varied,  the  re¬ 
sponse  function  and  burning  rate  were  both  increased  by  going  to  finer 
particle  size.  Bimodal  material  behaved  as  though  it  were  of  particle 
size  intermediate  to  its  ingredient  sizes  (with  possible  exception  at 
low  pressure  and  high  frequency  (Fig.  12(a)).  The  pressure  dependence 
of  y/e  was  greater  with  coarse  AP  than  fine  AP,  except  for  the  copper 
chromite  catalyzed  formulation  (Fig.  12(d)). 

The  enhanced  response  function  obtained  with  fine  oxidizer  is 
probably  due  to  an  increase  in  the  role  of  heat  release  steps  at  or  near 
the  burning  surface  due  to  high  interfacial  area  in  the  propellant  (Ref. 
14).  The  higher  pressure  sensitivity  with  coarse  oxidizer  is  probably 
related  to  the  pressure  dependence  of  the  AP  deflagration  in  the  200- 
800  psi  range  (Ref.  15),  and  the  correspondingly  important  role  of  the 
diffusion  flame.  The  absence  of  this  particle  size  effect  with  the 
formulation  containing  copper  chromite  is  presumably  related  to  the  as- 
yet-undetermined  mode  of  catalysis  by  copper  chromite,  which  not  only 
reduces  the  pressure  dependence  of  y/e,  but  also  has  a  major  depressing 
effect  on  y/e  at  all  pressures.  This  could  be  through  a  decoupling  of 
the  gat7  phase  reactions  and  the  surface  reactions,  but  the  mechanism 
remains  to  be  found. 

Effect  of  Ballistic  Modifiers 

Tests  were  run  on  three  different  families  of  propellants  involving 
ballistic  modifiers.  All  used  PBAN  binder.  One  family  contained  76%  of 
80y  AP  and  provides  a  comparison  of  behavior  with  1%  LiF,  1%  copper 
chromite,  and  no  modifier.  A  second  family  contained  76%  of  15y  AP  and 
again  provides  a  comparison  of  behavior  with  1%  LiF,  1%  copper  chromite, 
and  no  modifier.  The  third  family  contained  80%  bimodal  AP  (50-50  ratio 
of  micropulverized  and  as-received),  and  compared  behavior  with  0,  0.5, 
and  2.0%  copper  chromite. 
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bv  th^f/h2  CfSe ,°f  8°,U  AP|  the  y/e  values  were  fairly  well  correlated 
ZJk  /  abscifa  (Fig.  13(a)-(c) )  and  the  value  of  u/e  was  high  for 
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FIG.  13(a).  Effect  of  Ballistic  Modifiers  on  Response 
Function  (PBAN  Binder),  76%  80p  AP,  215  psia. 
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FIG.  13.  Effect  of  Ballistic  Modifiers 
Response  Function  (PBAN  Binder) . 
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One  formulation  was  prepared  as  a  modification  of  A-l  by  replace¬ 
ment  of  1%  of  binder  with  1%  of  AI2O3.  This  modification  produced  a 
small  decrease  of  y/e,  of  marginal  significance  (Fig.  14). 

Double-Base  Propellants 

Of  the  four  double-base  propellants  tested,  two  (N-4  and  N-5)  are 
of  roughly  the  same  energy  level  as  the  composites  tested,  while  two 
(JPN  and  X-14)  are  fairly  high  energy  and  high  burning  rate.  Two  (N-5 
and  X-14)  are  plateau  type  propellants  (i.e.,  they  exhibit  a  low  de¬ 
pendence  of  burning  rate  on  pressure  over  a  pressure  interval  in  the 
rocket  motor  operating  range) .  The  test  results  do  not  seem  to  reveal 
any  conspicuous  difference  in  response  functions  (see  Appendix  B) .  The 
values  of  y/e  are  somewhat  higher  than  for  most  of  the  composites,  but 
not  as  high  as  the  highest  given  with  composites.  The  values  of  y/e 
obtained  with  the  lower  energy  formulations  were  about  the  same  as  with 
the  higher  energy  formulations.  Of  course  the  high  energy  formulations 
are  much  more  unstable  (on  a  per  unit  burning  surface  area  basis)  be¬ 
cause  of  the  high  burning  rate  and  high  val^e  of  c  (Eq.  37).  The 
difference  between  the  plateau  and  nonplateau  propellants  was  not  con¬ 
spicuous  enough  to  belabor  in  the  absence  of  more  systematic  studies, 
although  the  plateau  propellants  had  somewhat  lower  values  of  y/e 
(particularly  in  the  case  of  N-5  versus  N-4  at  815  psia) .  It  is  notable 
that  the  variation  of  y/e  with  pressure  was  less  consistent  for  the 
double-base  propellants  than  with  the  composites. 

Effect  of  Pressure 

The  present  studies  have  provided  considerable  data  on  the  effect 
of  pressure  on  y/e,  as  evident  from  the  foregoing  discussion.  The  most 
complete  studies  of  pressure  dependence  are  provided  by  the  tests  on 
PBAN  formulations  A-13,  A-15,  A-17,  A-14,  A-16,  A-18  and  A-19  (Table  2, 
Fig.  12(c)  and  (d)  and  Fig.  15(a)-(c)  respectively).  As  noted  earlier, 
correlation  of  the  response  function  data  by  constant  pressure  curves 
plotted  against  f/r2  instead  of  f  yields  improved  similarity  among  the 
constant  pressure  curves,  as  suggested  by  current  one-dimensional  the¬ 
ories.  On  the  other  hand,  these  same  theories  do  not  generally  predict 
the  conspicuous  pressure  dependence  of  y/e.  The  nature  of  the  pressure 
dependence  was  the  same  for  all  the  composite  propellants  tested,  with 
high  response  function  at  low  pressure  and  declining  y/e  as  pressure 
increased.  The  degree  of  pressure  dependence  was  roughly  the  same  for 
all  formulations  except  A-14  and  A-19,  which  had  low  pressure  depen¬ 
dence  of  y/e.  These  two  formulations  were  without  ballistic  modifiers 
and  were  effectively  fine  AP  (and  high  burning  rate,  high  y/e).  As 
suggested  earlier,  the  low  pressure  sensitivity  may  be  due  to  the  rela¬ 
tively  high  importance  of  exothermic  surface  reactions  related  to  the 
high  interfacial  area  between  binder  and  oxidizer. 
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Effect  of  Additives  on  Attenuation 

The  effect  of  condensed  phase  products  has  been  discussed  in  many 
references,  the  most  appropriate  for  the  present  purpose  being  Ref.  16 
(which  shows  calculated  «2  versus  frequency  for  several  particle  sizes) . 

In  the  present  work,  the  condensed  phase  material  presumably  consists  of 
copper  chromite,  Al,  Al-^Oj  and  carbon.  Particle  sizes  were  not  measured. 
Copper  chromite  is  roughly  1  to  2  y  in  diameter  when  added  to  the  pro¬ 
pellant.  Aluminum  size  (burning  droplets)  is  presumably  of  the  same 
order  as  the  ingredient  material.  AI2O3  reaction  products  are  presumably 
either  in  the  <  2  p  range  (smoke  product) ,  or  of  a  size  related  to  the 
original  Al  or  agglomerates  (products  of  surface  reaction  on  Al  droplets)  . 
Carbon  is  probably  submicron  in  size. 

Looking  at  the  results  in  Fig.  9,  the  02* s  for  A-13  and  A-14  are 
lowest,  as  would  be  expected  since  no  additives  were  present.  The  effect 
of  addition  of  copper  chromite  can  be  estimated  by  comparison  of  results 
for  A-13  or  A-14  with  A-l  (1%  copper  chromite)  ,  which  led  to  a  40%  in¬ 
crease  in  02*  This  result  is  important  because  all  of  the  other  curves 
except  one  showing  effect  of  additives  include  the  presence  of  1  to  2  y 
copper  chromite. 

The  largest  attenuations  were  obtained  with  A-2  propellant,  con¬ 
taining  1%  5  y  Al  and  2%  copper  chromite.  Experience  with  aluminum  com¬ 
bustion  suggests  that  5  y  Al  would  burn  quickly  to  'v  1  y  AI2O3  smoke, 
which  is  a  relatively  good  size  for  attenuation  (see  Fig.  3  of  Ref.  16). 
The  value  of  02  *s  somewhat  low  but  in  qualitative  agreement  with  the 
theory.  The  result  with  0.5%  5  y  Al  (A-5)  lies  halfway  between  the  re¬ 
sults  for  1.0%  Al  and  0.0%  Al  (A-l  and  A-2).  A-20  provides  an  indica¬ 

tion  of  the  effect  of  adding  1%  of  20  y  Al  to  a  formulation  with  no 
copper  chromite  (e.g.,  A-14).  The  increase  in  02  is  not  very  great — 
possibly  because  20  y  Al  is  too  coarse  for  effective  attenuation  and  may 
not  react  fast  enough  to  convert  to  finer  AI2O3.  Addition  of  1%  20  y 
Al  to  a  formulation  with  2%  copper  chromite  had  a  more  conspicuous 
effect  (compare  A-7  and  A-14),  although  less  effect  than  5  y  Al.  Addi¬ 
tion  of  1%  2  y  AI2O3  (A-ll)  and  of  0.4%  0.02  y  Al  had  negligible  effect 
on  ci2  (compare  with  A-l) .  This  result  is  surprising  in  the  case  of  the 
AI2O3,  since  the  particle  size  was  near  optimum  for  damping:  it  is 
possible  that  the  oxide  was  agglomerated  to  an  ineffective  size,  possibly 
during  propellant  processing.  No  direct  information  on  this  point  is 
available. 

Stability  Limits 

In  the  present  work,  considerable  effort  was  devoted  to  testing 
near  the  stability  limits  of  the  combustor — primarily  to  extend  the 
range  of  y/e  data.  In  the  process,  considerable  information  regarding 
the  stability  limit  of  this  type  of  T-burner  was  developed,  and  curve¬ 
fitting  procedures  were  contrived  to  maximize  the  generalization  of 
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these  results.  The  ultimate  application  of  these  results  could  be 
either  (a)  as  an  empirical  guide  to  choice  of  test  variables  in  future 
work,  or  (b)  development  or  evaluation  of  loss  and  scaling  laws  for 
T-burners.  The  former  of  these  applications  can  be  made  from  the  re¬ 
sults  as  they  are  presented  (Table  5  and  Appendix  A) — at  least  for 
similar  propellants  and  T-burners.  The  second  application  is  beyond 
the  scope  of  the  present  study,  but  the  data  may  be  a  stimulus  to  such 
studies. 

It  should  be  clear  that  the  curve  fitting  procedure  that  was 
chosen  was  a  matter  of  expediency  in  the  face  of  limited  data.  The 
resulting  stability  limit  curves  differ  markedly  as  a  function  of  pro¬ 
pellant,  and  consequently  can  be  used  as  qualitative  indicators  of  the 
stability  "personality"  of  the  propellant.  In  practice,  of  course,  one 
would  be  mc’-e  concerned  with  how  high  the  response  function  becomes 
inside  the  stability  limit,  but  the  limit  is  an  indicator  of  the  range 
of  frequency  and  pressure  over  which  instability  might  occur. 


0  4  0  80  120  160  200 

f /re 


(a)  76%  80y  AP. 

FIG.  15.  Effect  of  Pressure  on  Response  Function 
(PBAN  Binder,  No  Ballistic  Modifier). 


Appendix  A 


NwC  TP  4393 


m 

CN 

o 

rH 

rH 

cn 

m 

o 

• 

CN 

• 

vO 

u-i 

m 

M  Pb  o  US 


fr-5 

i  ^ 

O 

« 

00 

rH 

•l 

m 

/~n  e-e 

c 

•*  rH 

S3  cn 

o 

•H 

- 

PQ 

^  rv 

4J 

•H 

W  3. 
ph  m 

AP 

cc 

CO 

O 

■  • 

•  « 

a 

e 

•  • 

•  • 

o 

g 

Hi  4J 

CJ 

•  p 

a)  co 

w  c 

N  tN 

ai  -h 

•H  »H 

•a  b 

TJ  CO 

C  3 

iH  iJ 

•H  fH 

X  CO 

PQ  < 

O  U 

lo  <r 

CN  O 

rH 

•h  cn 

m  o 

• 

CN  • 

vo  m 

m 

CN 

CO 

00 

4J 

cn 

CO 

•  • 

•  • 

• 

Q 

CO  • 

•  • 

• 

O 

•H 

•H 

CO  • 

•  • 

• 

U 

a 

CO 

• 

•  • 

• 

•H 

*  o 

o 

rH 

a;  <u 

•  V 

• 

rH 

Hi  CO 

CO 

CO 

P  \ 

•  \ 

• 

PQ 

CO  c 

a 

CO  *H 

•  «H 

• 

a) 

Hi 

M  •'* 

•* 

PH  U 

Ph  a 

<,-»••>  D 


B  b  ti 
•  9  <U  Cfl 
[j  6  N  >1 

fl]  *H  Tl  H 

TJ  g  TJ  « 
B  3  iH  u 
rj  x  « 
«  <  o  u 


yst.  .  cc,  2%  c,  in/sec.  .  .  .38,500  Catalyst.  .  cc,  2%  c,  in/sec.  .  .  .38,  500 

Kr . 25.1  Kr .  36.7 


ALUMINIZED  COMPOSITE  PROPELLANTS 


Appendix  B 


EXTENDED  PRESSURE  COVERAGE 


Propellant  (17)  A-13 


Composition 

Ballistic  Data 

Binder.  .  .  .  PBAN,  24% 
Aluminum.  .  .None 

Oxidizer.  .  .  AP(80y) ,  76% 
Catalyst.  .  .  None 

Pressure , 
r,  in/sec 
P/r  .  .  . 
c,  in/sec 
Kr.  .  .  . 

psia.  .  .  215 

.  •  16 

. 1344 

. 35,400 

. 64.1 

415 

.21 

1976 

36,700 

91.2 

815 

.28 

2911 

37,600 

130.0 

63 


Response  Function  Results  and  Neu 
tral  Boundary  for  A-13  Propellant 


M 


80 

a2 


Propellant  (18)  A-15 


Composition 

Ballistic  Data 

Binder.  .  .  .  PBAN,  23% 

Pressure, 

psia.  .  .  215 

415 

815 

Aluminum.  .  .  None 

r,  in/sec 

•  •  •  •  •  • 26 

.35 

.47 

Oxidizer.  .  .  AP(80p),  76% 

P/r  .  .  . 

.  827 

1186 

1734 

Catalyst.  .  .  cc,  1% 

c,  in/sec 

. 36,500 

37,200 

37,700 

Kr .  .  .  . 

. 37.9 

54.1 

77.6 

Response  Function  Results  and  Neu¬ 
tral  Boundary  for  A-14  Propellant 
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Composition 


Binder.  . 
Aluminum. 
Oxidizer. 
Catalyst. 


PBAN,  23% 
None 

AP(15u),  76% 
cc,  1% 


Ballistic  Data 


r,  in/sec 
P/r  .  .  . 


c,  in/sec 


Kr. 


psia.  .  .  215 

415 

. 51 

.64 

.  422 

648 

. 37,000 

38,200 

38 

. 19.3 

28.7 
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Propellant  (24)  A-19 


Composition 

Binder.  .  .  .  PBAN,  24% 

Aluminum.  .  .  None 

Oxidizer.  .  .  AP(80y),  76% 

Catalyst.  .  .  LF,  1% 


Ballistic  Data 


Pressure,  psia.  . 

215 

815 

r,  in/sec  .  .  .  . 

.21 

.34 

P/r . 

.  1024 

2397 

c,  in/sec  .... 

.35,900 

37,600 

Kr . 

.  48.2 

108.2 

1615 

.44 

3670 

38,100 

163.1 
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Response  Function  Results  and  Neu 
tral  Boundary  for  A-19  Propellant 


Response  Function  Results  and  Neu 
tral  Boundary  for  T-35  Propellant 


f  (kcps) 


Propellant  (34)  JPN 


Ballistic  Data 


in/sec 


Response  Function  Results  and  Neu¬ 
tral  Boundary  for  X-14  Propellant. 
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Appendix  C 

CHARACTERIZATION  OF  THE  STABILITY  LIMIT  AND 
[3(u/e)/3f]  ,  [3(y/e)/3p]f  =  0  CURVES 

Experience  with  the  1.5-inch  diameter  T-burner  has  shown  that,  at 
any  given  mean  pressure,  there  is  a  limited  frequency  range  in  which 
04  >  0,  and  at  any  given  frequency  (burner  length)  there  is  a  limited 
pressure  range  in  which  04  >  0.  While  the  cost  and  accuracy  of  deter¬ 
mination  of  the  a  =  0  stability  limit  has  previously  discouraged  de¬ 
termination  of  the  limit,  the  present  investigation  has  explored  this 
subject  as  thoroughly  as  the  need  seems  to  warrant.  This  investigation 
included  an  analytical  representation  of  the  stability  limit  designed  to 
extract  as  much  information  as  possible  from  the  rather  limited  data 
available. 

The  stability  limit  is  characterized  by  upper  and  lower  pressure 
limits  Pc  and  PG  beyond  which  the  burner  is  stable  for  all  lengths 
(frequencies)  (Fig.  7).  Likewise  there  are  upper  and  lower  frequencies 
fc  and  f0,  beyond  which  the  burners  are  stable  at  all  pressures.  Join¬ 
ing  the  points  Pc  and  P0  is  a  curve  (designated  fp  in  Fig.  7)  indicating 
the  frequencies  for  maximum  response  function,  at  each  pressure, 
[3(p/e)/3f]p  =  0;  and  a  corresponding  curve  (designated  ff)  joining 
points  fn  and  fm  is  the  locus  of  pressures  for  maximum  p/e  for  each  fre¬ 
quency,  [3(y/e)/8p]f  =  0.  The  intersection  of  these  curves,  in  the  un¬ 
stable  region  (point  designated  P*,  f*)  is  the  condition  for  maximum 
response  function. 

A  study  of  stability  limit  data  obtained  from  test  series  at  various 
mean  pressures  showed  that  the  upper  frequency,  fu,  for  x  •=  0  could  be 
correlated  with  pressure,  by  a  linear  relation  between  fu/r2  and  Pu. 
Accordingly,  this  functional  form  was  adopted  for  correlation  of  all 
upper  frequency  limit  data.  The  detailed  form  was  modified  for  conve¬ 
nience  to  represent  the  upper  and  lower  absolute  pressure  limits,  Pc  and 
PQ,  and  the  upper  frequency  limit,  fm,  as  singular  points: 
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The  forms  showing  Pm  and  fm  as  singular  points  are  determined  by  solving 
Eq.  41  for  P,  setting  SP/3f  =  °°,  and  using  the  result  to  obtain  Pm  from 
the  equation  for  P.  The  result,  using  also  the  burning  rate  relation 
r  =  CPn,  is 


2  2 
f  r  P  -f  r  P 

o  C  C  CO  o 


/  _  V  ft  r  P  -f  r 
|  2n  \  /  o  c  c  c 

"  \ l+2n  j  2  . 

'  '  f  r  -f  r 


o  c  co 


which  is  approximated  by  the  simpler  form 


P  =  ( -  -2n-  |  p 
m  \  l+2n  I  c 


2  2  2  1 
for  the  combinations  of  fcrQ  PQ  with  fQrc  Pc  and  of  fcrG  with  f0nc* 

encountered  in  this  work.  Corresponding  forms  for  the  frequency  at 

maximum  are 


r  /  f  r  \  <  P  . 

f  _  _HL  f  -2 - 

m  2n  1  2  2 )  \  P  -P  / 

-v*  '  r*  r\  ' 


m  l+2n  \  r 


A  study  of  the  experimental  data  for  the  low  frequency  limit  data 
lead  to  correlation  by  the  function 


/p-p  r 

Mr) 


Most  of  the  experimental  f^  values  seem  to  indicate  an  exponent  for  the 
pressure  terms  somewhat  smaller  than  3,  but  taking  into  account  the  un¬ 
certainty  due  to  data  scattering,  the  exponent  3  provides  a  more  con¬ 
servative  estimate.  The  terminal  point  (Pn,  fn) ,  where  fn  is  the  lowest 
frequency  sustaining  unstable  burning  and  Pn  the  corresponding  pressure, 
is  calculated  by  applying  the  condition  for  minimum,  df^/dP  =  0,  to 
Eq.  46.  The  equation  for  Pn  becomes 

P  [(f  f  -  f  ]  +  Pn  [f  -  (f  f  )**] 
n - <47) 

c  o 
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and  the  approximate  solution  for  fn  is  given  by 

k  3 


f  " 

C  0  c 


f  -f 
c  o 


+  f 


k  3 

f  f  r  -  f 

O  C  0 


f  -f 
c  o 


(48) 


The  locus  lines  [3(p/e)/3f]g  *  0  and  [3(p/e)/3p]£  =  0  were  assumed 
(after  examination  of  available  data)  to  be  of  the  form 

fp-fo(t)(?V)  +  £c  (f-)(£i-) 


(49) 


and 


n  m  n  '  m '  m  n 


(50) 


where  the  functional  form  has  been  chosen  so  that  the  curves  go  through 
the  points  (fc,Pc)  and  (f0,P0)  in  one  case  and  (fn,Pn)  and  (Pm,Pm)  in 
the  other  case.  Sample  curves  are  shown  in  Fig.  6. 

The  intersection  of  the  [3(p/e)/3f]  -  0  and  [3(p/e)/3P]  ■  0 

curves  is  the  point  of  maximum  response  function,  denoted  here  by 
(f*  P*).  This  intersection  could  be  calculated  by  graphical  means,  or 
estimates  could  be  made  with  the  approximate  relations 


f* 


1+n  r*  f 

2  r  o 
o 
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MAXIMUM  ACOUSTIC  PRESSURE  (AP)m,  PSIA,  AND  (AP)m/P  RATIOS 
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